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ABSTRACT 

To  support  the  design  and  analysis  of  passive  infrared  atmospheric  sounding  instruments,  an 
analytical  error  analysis  technique  has  been  implemented.  This  technique  is  based  on  a  linear  approxi¬ 
mation  to  the  radiative  transfer  equation  and  uses  a  minimum  variance  estimation  approach  to  atmo¬ 
spheric  profile  retrieval.  As  a  result,  once  the  proper  matrices  have  been  constructed  an  estimate  of  the 
retrieval  error  is  computed  through  a  single  linear  matrix  equation.  The  solution  vector  is  written  with 
temperature  and  water  vapor  explicitly  defined,  thus  producing  the  error  vector  for  both  simultaneously. 
Examples  showing  typical  results  are  presented  for  a  high  spectral  resolution  sounder  as  well  as  for  a 
lower  resolution,  filter-wheel  type  of  instrument. 


j  Accesion  For  "j 

NTIS 

CRA&I 

DTIC 

r  A  3  .’"i 

IT;  1  ■ !  V.  -J  r 

j  J'-itifiCcit.Ofi 

By_ 

Dhd.  lb 

litiO.I  / 

Dist 

MV3;l  tj'.fj/or 

S;.-.rC'al 

/)'/ 

i 

1  . 

iii 


ACKNOWLEDGMENT 


Henry  E.  Fleming  is  acknowledged  for  his  contributions  to  the  science  of  satellite  atmospheric 
sounding,  including  the  publication  of  the  retrieval  error  algorithm  on  which  this  work  is  based. 


TABLE  OF  CONTENTS 


Abstract  iii 

Acknowledgment  v 

List  of  Illustrations  ix 

List  of  Tables  xi 

1.  INTRODUCTION  1 

1 . 1  Motivation  1 

1.2  Context  and  Limitations  1 

1.3  Report  Overview  2 

2.  THEORY  AND  IMPLEMENTATION  3 

2.1  Radiative  Transfer  Equation  3 

2.2  Numerical  Implementation  5 

2.3  Estimation  of  Retrieval  Operator  and  Error  Covariance  12 

3.  A  PRIORI  TEMPERATURE  AND  WATER  VAPOR  STATISTICS  15 

4.  EXAMPLE  RETRIEVAL  ERROR  ANALYSES  21 

4. 1  Introduction  2 1 

4.2  AIRS  21 

4.3  GOES-I  Sounder  24 

5.  SUMMARY  27 

APPENDIX  A.  FASCOD3  INPUT  29 

LIST  OF  SYMBOLS  39 


REFERENCES 


41 


LIST  OF  ILLUSTRATIONS 


Figure 

No.  Page 

1  Triangular  spectral  response  array  for  case  with  0.05-cm'’  bandwidth.  1 1 

2  Mean  temperature  profiles  for  January  (top)  and  August  (bottom).  17 

3  Mean  water  vapor  profiles  for  January  (top)  and  August  (bottom).  18 

4  Standard  deviation  of  temperature  for  January  (top)  and  August  (bottom).  19 

5  Standard  deviation  of  water  vapor  for  January  (top)  and  August  (bottom).  20 

6  AIRS  retrieval  error  for  temperature  (top)  and  water  vapor  (bottom).  23 

7  GOES-I  sounder  retrieval  error  for  temperature  (top)  and 

water  vapor  (bottom).  26 


IX 


LIST  OF  TABLES 


Table 


No. 

Page 

1 

Pressure  Levels  Used  in  Analysis 

15 

2 

Month  and  Latitude  Statistics 

16 

3 

Model  AIRS  Instrument  Specifications 

22 

4 

Model  AIRS  Spectral  Regions 

22 

5 

Model  GOES-I  Sounder  Specifications 

24 

6 

Model  GOES-I  Sounder  Channels 

25 

XI 


1.  INTRODUCTION 


1.1  MOTIVATION 

The  science  of  atmospheric  temperature  and  water  vapor  profile  retrieval  from  space-based  pas¬ 
sive  infrared  and  microwave  sensors  has  been  developing  since  the  i^60s  [1,2]  and  is  routinely  applied 
to  satellite  measurements  for  use  in  weather-forecasting  [3]  and  climate-monitoring  [4]  applications. 
However,  the  accuracy  and  resolution  of  the  retrieved  profiles  from  current  generation  instruments  are 
such  that  the  use  of  their  data  leads  to  mixed  results  [5]  and  falls  short  of  what  is  desired  from  future 
instruments  [6].  These  next-generation  sounders  must  be  designed  to  collect  data  that  will  result  in 
significant  improvements  in  the  quality  of  the  retrieved  profiles. 

The  analysis  of  the  expected  quality  of  retrievals  from  an  instrument  design  is  in  general  a 
difficult  task.  Whereas  the  raw  data  quality  from  an  instrument  (e.g.,  signal-to-noise  ratio,  spatial,  and 
spectral  resolution)  can  be  objectively  calculated  given  a  specific  design,  the  retrieval  process  depends 
on  algorithms  and  sources  of  data  external  to  the  instrument.  Retrieval  algorithms  vary  from  those 
specifically  based  on  a  statistical  regression  between  collocated  in-situ  measurements  and  satellite  sound¬ 
ings  to  those  which  use  a  radiative  transfer  model  along  with  linear  estimation  theory  and  iterate  from 
a  first  guess.  The  quality  of  these  retrievals  depend  in  a  large  part  on  the  a  priori  knowledge  of  the  state 
of  the  observed  atmosphere.  These  variabilities,  along  with  the  fact  that  data  collected  are  used  for  other 
applications  such  as  trace  gas  monitoring  and  the  calculation  of  cloud  heights  and  amounts,  challenge  the 
characterization  of  “performance”  for  atmospheric  sounding  instruments. 

Given  these  challenges,  however,  it  still  remains  an  important  task  to  evaluate  the  expected 
performance  of  a  sounding  instrument  both  to  justify  the  instrument  design  requirements  as  well  as  to 
judge  the  impact  of  design  modifications  that  may  become  necessary  during  instrument  development.  The 
most  accurate  way  of  predicting  on-orbit  performance  is  a  full  simulation  of  a  variety  of  atmospheric 
situations  and  the  application  of  the  operational  retrieval  algorithm  to  the  resulting  data.  This  prediction 
would  have  to  be  done  for  every  point  in  the  design  trade  space  and  would  not  be  desirable  given  finite 
resources. 

Alternatively,  a  simplified  approach  to  retrieval  error  analysis  based  on  linear  estimation  theory 
can  be  implemented  as  described  in  this  report.  This  algorithm,  based  on  previous  work  [7],  uses  a  linear 
retrieval  model  to  estimate  the  variance  of  the  error  in  the  retrieved  temperature  and  water  vapor  profiles. 
Since  the  algorithm  can  be  described  with  an  analytical  equation,  parameter  sensitivity  studies  can  be  run 
quickly  and  aid  in  the  development  of  a  sounding  instrument. 

1.2  CONTEXT  AND  LIMITATIONS 

The  application  and  results  of  the  methodology  documented  in  this  report  should  be  considered 
in  light  of  its  limitations.  First,  it  is  not  a  full  retrieval  analysis.  While  the  basic  retrieval  operator  derived 
could  be  applied  to  simulated  or  measured  data,  this  report  focuses  on  its  use  in  the  estimation  of  expected 
retrieval  error.  Thus,  the  algorithm  does  not  require  sounding  data  but  rather  a  description  of  the  instru- 
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merit  including  the  number,  bandwidth,  noise  level,  and  spectral  location  of  the  channels.  No  spatial 
effects  are  considered  including  off  axis  sounding  or  the  impact  of  clouds.  Some  of  these  limitations 
could  be  remedied  with  extensions  to  the  work  presented. 

Fundamentally,  the  error  analysis  algorithm  involves  a  linear  expansion  around  the  known  solu¬ 
tion  to  the  nonlinear  retrieval  problem.  This  formulation  has  the  impact  of  underestimating  the  retrieval 
error  in  most  situations.  The  results  should  be  considered  to  be  a  lower  bound  on  the  error.  This  char¬ 
acteristic  of  the  results  leads  naturally  to  the  question;  if  the  error  levels  predicted  are  optimistic  compared 
to  full  retrievals,  what  can  be  said  about  the  sensitivity  to  instrument  parameters?  This  issue  has  been 
considered  and  while  no  guarantee  can  be  provided,  the  impact  on  retrieval  error  of  instrument  parameter 
variations  can  be  evaluated  in  a  relative  sense  with  the  general  trends  expected  to  hold. 

In  addition  to  the  above  limitations  of  this  approach,  it  should  be  noted  that  only  temperature  and 
water  vapor  are  explicitly  considered  in  the  retrieval  analysis.  The  performance  of  the  instrument  design 
in  other  sounding  applications,  such  as  trace  gas  monitoring,  has  not  been  included. 

The  most  appropriate  application  of  this  methodology  is  in  the  relative  comparison  of  various 
design  options  during  a  top-level  system  analysis  study.  The  ease  of  implementation  and  use  as  well  as 
speed  in  computation  allow  a  large  number  of  design  options  to  be  considered.  After  a  design  has  been 
constrained,  final  design  choices  can  be  made  using  performance  predictions  based  on  full  simulations 
and  retrievals. 

1.3  REPORT  OVERVIEW 

The  main  goal  of  this  report  is  to  document  in  detail  the  algorithm  used  to  estimate  the  error  in 
retrieved  atmospheric  temperature  and  water  vapor  profiles  given  an  instrument  description.  This  algo¬ 
rithm  has  been  implemented  at  Lincoln  Laboratory  and  has  formed  the  basis  for  several  sounding  system 
studies.  Section  2  provides  the  details  behind  the  theory  and  implementation  of  the  technique.  Section  3 
presents  mean  vectors  and  covariance  matrices  of  atmospheric  temperature  and  water  vapor  profiles 
obtained  from  NOAA.  These  statistics  constitute  the  a  priori  information  used  in  the  error  analysis 
technique.  Section  4  presents  some  example  results  of  the  technique  applied  to  two  different  sensors. 
Section  5  concludes  the  report  with  a  summary. 
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2.  THEORY  AND  IMPLEMENTATION 

2.1  RADIATIVE  TRANSFER  EQUATION 

The  fundamental  equation  for  atmospheric  sounding  from  a  satellite  is  the  radiative  transfer 
equation  (RTE)  shown  in  Equation  (2.1). 

1 

R=B  X  +  [adz  ,  (2.1) 

r 

S 

where  s  indicates  surface  conditions  and 

R  =  upwelling  radiance  seen  by  satellite  ir  a  particular  spectral  channel  in  mW/(m^-Sr;, 

B  -  radiance  emitted  by  a  given  layer  of  the  atmosphere  in  a  particular  channel,  and 

T  =  transmittance  from  en  layer  of  atmosphere  to  top  of  atmosphere  for  a  particular  channel. 

The  upwelling  radiance  R  is  in  genera)  a  function  of  temperature,  pressure,  the  concentrations  of 
various  atmospheric  constituents  as  well  as  the  central  wavenumber  and  spectral  response  of  the  channel. 
In  this  analysis  the  concentration  of  all  atmospheric  constituents  except  water  vapor  is  assumed  to  be 
known  and  constant.  Also,  sensor  noise  will  be  neglected  for  now  but  will  be  added  in  later. 

Considering  variability  of  atmospheric  temperature  and  water  vapor  profiles,  a  mean  state  can  be 
defined  as  in  Equation  (2.2). 


R  =  B  X  +  {BdX  (2.2) 

5  S  J 

Assuming  that  the  radiance  and  the  transmittance  are  uncorrelated  (not  entirely  true,  but  a  prac¬ 
tical  assumption).  Equation  (2.2)  can  be  rewritten  as  follows. 

1 

R  =  WT+{bJx  (2.3) 

S  S  J 

r 

s 

Now,  defining  the  difference  quantity  to  be  the  mean  subtracted  from  the  true  value  (i.e., 
Ax  =  x-x)  subtracting  Equation  (2.3)  from  Equation  (2.1)  results  in  Equation  (2.4). 
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I  1 

AR  =B  T  -B  r  +  \Bdr-  { Bdr 

S  S  S  S  J  J 

T  T 

s  s 

=  B  T  -W T 

5  S  S  S  S  S  S  S 


+  j[Bdr- Bdr  + Bdr-Bdr) 


I  I 

=  AB  T  +  B  At  +  [ABdx+  [BAdr 

s  s  s  s  }  J 


(2.4) 


By  integrating  the  last  term  of  Equation  (2.4)  by  parts,  the  following  equation  results. 

1  140- 

^BAdx  =BAx  -^\x~dX 


h  1 


dX 


40 

=  0-F Ax  -  f 

s  s  J 


1 


dB 
Ax~dX 
dX 


(2.5) 


Here,  an  auxiliary  variable  X  is  used  which  represents  the  layers  of  the  atmosphere.  The  pressure 
levels  at  the  bottom  of  each  of  the  40  layers  used  in  the  analysis  are  defined  in  Section  3.  Substituting 
Equation  (2.5)  into  Equation  (2.4)  results  in  Equation  (2.6). 

40  HR 

AR  =  AB  X  +  \AB—dX-  \Ax~dX  (2.6) 

^  ^  J  dX  J  dX 

I  I 

The  radiance  of  a  layer  depends  on  the  layer’s  temperature  through  the  nonlinear  Planck  function, 

but  one  can  approximate  the  radiance  difference  factor  AB  by  the  first  term  of  a  Taylor’s  expansion  about 
the  mean.  This  expansion  is  shown  in  Equation  (2.7)  where  At  is  the  difference  between  the  true  tem¬ 
perature  of  that  layer  and  its  mean  temperature. 


dB  ^ 

AB  =  —  At 
dt 

Substituting  Equation  (2.7)  for  the  AB  terms  in  Equation  (2.6)  yields  Equation  (2.8). 


(2.7) 


dB 


40 


.  r.  ^  dB  dx  r  ^  dtf 

AR  = — ^At  X  +  Ar - dX-  At — dX 

dt  ^  ^  J  dt  dX  J  dX 


40  — 

dB 


(2.8) 
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So  far,  only  the  effect  of  temperature  has  been  considered  explicitly.  The  effect  of  water  vapor 
on  the  upwelling  radiance  is  included  by  expanding  the  At  in  the  last  term  of  Equation  (2.8)  about  the 
mean  value  of  precipilable  water  u  measured  in  centimeters. 


At  =  — 


du  l“=“ 


Au 


) 

du 

dr. 


Au 


D 

T  - —+  T  - ^ 

^  du  ^  du 


Au 


0+  T 


dr. 


w 


®  du 
rw  du 


Au 


Au 


(2.9) 


Here,  T  =  t^Tj^,  and 

=  transmittance  of  dry  atmosphere  with  no  precipitable  water,  and 
r^  =  transmittance  of  wet  atmosphere  with  only  precipitable  water. 

Equation  (2.9)  can  then  be  substituted  into  Equation  (2.8)  to  result  in  Equation  (2.10),  which 
explicitly  shows  the  upwelling  radiance  as  a  function  of  atmospheric  temperature  and  precipitable  water. 


dB 

AR  =  — ^  T  At 
dt  ^  " 


+ 


40  — 

f  dB  dr 

Ar - 

J  dt  dX 

1 


dB 
du  dX 


AudX 


(2.10) 


Note  that  precipitable  water  is  defined  to  be  0  above  the  1 5th  layer  in  the  atmosphere. 


2.2  NUMERICAL  IMPLEMENTATION 
2.2.1  Numerical  Quadrature 

The  integrals  in  Equation  (2.10)  can  be  evaluated  numerically  by  a  variety  of  quadrature  rules. 
In  this  analysis,  the  following  equation  is  derived  for  two  arbitrary  functions  fiX)  and  g(X). 
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40  40  y+>  . 

\f{X)g{x\~dx  =Y  f  f{X)giX)^dx 

J  dX  “  ^  dx 

1  y=l  j 

40  y+> 

=  y7i”  f  — 

}  J  dx 
y=i  > 


Y  -Q  ■(/>»  ^7-t-i 

y-.  2 


T.  ,  -  T. 

y+1  y 


40  41 

~  2  ^  ■^y  ^y  1  ^y+'  ~  ^y ]  ^  2  ^  ^-'l 1 

y=i  y=2 

•'1*1  9  -^■'y^y  2 

L  ^  J  J=2  L  ^  J 

^/41^4l[-^^^  • 


(2.11) 


Here,  note  that  the  pressure  levels  are  defined  as  in  Table  1.  Also,  T^,  =  1  and/4,  ^41  “/40  ^40 
since  the  41st  layer  is  above  the  atmosphere.  Thus,  using  Equation  (2.11),  Equation  (2.10)  can  be 
rewritten  in  a  discrete  form  as  follows. 


40  15 

AR.  =  b.At  +  Y  b.  A/ .  -  y  c..  A« . 

I  SI  S  ^  IJ  J  ^  tj  J 
7=1  7=1 

The  index  i  refers  to  the  spectral  channel.  The  coefficients  are  defined  as  follows. 


(2.12) 


6 


b  . 

SI 


f  dB  ^ 


V  y 


T  . 

SI 


r 

[dt, 

il 

2 

(  dB^ 

\j+l  ^ij-l 

\  > 

ij 

2 

f  dB" 

'  2  —  T  —  T 
i40  t39 

i40 

2 

(  X  dx^^ 

1 

_ 1 

it 

L  2  J 

(  r 

B.  .  ,-B.  .  , 

(.7+1  i,j-\ 

y 

2 

for  j  =  1 


for  2  <  y  <  39 


for  7  =  40 


for  7  =  1 


for  7  =  2  to  15 


The  computation  of  dx^/du  is  shown  later  in  Section  2.2.4. 

2.22  Conversion  From  Precipitable  Water  u  to  Miidng  Ratio  q 

The  water  vapor  statistics  that  are  used  in  this  analysis  are  given  in  terms  of  water  vapor  mixing 
ratio  q  (g/kg)  rather  than  precipitable  water  u.  Thus,  a  change  of  variable  for  the  last  term  in  Equa¬ 
tion  (2.12)  is  necessary.  The  amount  of  precipitable  water  u  in  the  atmosphere  above  a  pressure  level  p 
is  given  by  Equation  (2.13). 

j  P 

“{p)  =  -l<}{p')dp'  ■  (2.13) 

^0 

Here,  g  =  980  g/(cm-sec^)  is  the  gravity  constant.  By  taking  delta  quantities,  reversing  the  limits 
of  integration,  and  applying  the  quadrature  result  of  Equation  (2.11),  the  following  equation  results. 
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15  ,  P*+l 

=  -%-  \^{p')<ip' 
k=J^  P, 

^ k=j^  ^ 


k+\ 


-Pk) 


J_ 

2« 


16  15 

S  ^k  (Pk-\  -Pkh'L^kiPk-  Pk^l ) 

k=j+\  k=j 


“  i*’;  ‘  '’jti  )■" i ^»  ('’*-!  ■  '’<*1 ) 


2g  ^ 

^  k=j+\ 


15 


2j?  ^ 

k=j 


where 


^k=Pk-Pk^i  for  k=j 

=  Pk-i-PM  fork>j 

Using  Equation  (2.14),  the  last  term  of  Equation  (2.12)  can  be  replaced  as  follows. 


15 

I 

y=i 


15 


j=l 


15 


*=1 


15 


28 


T^ij^k 

Lv=i 


Let 


^ik  2e  ^ 
«>=1 


(2.14) 


(2.15) 


(2.16) 
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where  is  as  defined  above. 

Then,  Equation  (2.12)  can  now  be  written  as  in  Equation  (2.17). 

40  15 

AR.=b  At  +y  b..At.-y  d  Ag.  (2.17) 

I  SI  S  ^  IJ  J  ^  IJ 

j=l  j=l 

This  shows  the  explicit  dependence  on  temperature  and  water  vapor  in  the  forward  radiative  transfer 
equation. 

2.23  Matrix  Solution  to  RTE 

Equation  (2.17)  can  be  written  in  matrix  form  as  Equation  (2.18). 

b  =  Av  (2.18) 

Here,  b  is  the  radiance  difference  vector,  A  is  the  radiative  transfer  matrix,  and  v  is  the  solution  difference 
vector  as  below.  N  is  the  number  of  spectral  channels. 


V  = 


^2 

^^40 

At 

s 

Aq^ 


15 


b,  . 

•  •  b. 

b  ,  1 

•  •  d, 

1.1 

• 

• 

• 

1,40 

• 

sA  1 

• 

• 

• 

1.1 

• 

• 

• 

1,15 

• 

• 

• 

• 

• 

b 

•  •  b 

• 

• 

b  1 

• 

• 

d 

• 

• 

•  •  d  , 

_NA 

^N,40 

s.N  1 

NA 

NA5] 
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This  represents  a  linearized  version  of  the  forward  RTE.  Section  2.3  discusses  solutions  to  this 
equation  given  measured  radiances  as  well  as  the  computation  of  the  covariance  of  the  retrieval  estimates. 

2.2.4  Calculation  of  Channel  Radiances  and  Transmittances 


The  channel  radiances  and  transmittances  used  above  include  the  effect  of  the  instrument  channel 
spectral  response.  This  section  details  their  computation. 

Channel  Radiances.  The  spectral  radiance  emitted  by  a  layer  with  temperature  T  (Kelvin)  at 
wavenumber  v  is  assumed  to  have  a  blackbody  spectrum  with  unit  emissivity  and  thus  is  computed  with 
the  Planck  function. 


where 


''  m^-Sr-cm~’ 

£  ^  -1 


(2.19) 


C,  =  1.191062x10-5 


K 


mW 


m^-Sr-^cm 


C,  =  1.438786  , 

^  cm” 


and 


are  the  appropriate  constants. 

The  derivative  of  the  Planck  spectral  radiance  function  with  respect  to  temperature  is  shown  in 
Equation  (2.20). 


dB 

_ V 

dT 


C  CjV 


T  _1 


mW 

m^  -  Sr -cm”*  -K 


V  J 


(2.20) 


The  bandwidth  of  a  channel  is  assumed  to  be  small  enough  that  the  Planck  function  and  its 
derivative  can  be  assumed  to  be  constant;  thus,  the  in-channel  quantities  of  the  spectral  radiance  and  its 
derivative  are  found  by  multiplying  them  by  the  bandwidth  of  the  channel.  The  radiance  factors  used  in 
Equation  (2. 1 2)  are  computed  for  channel  /  with  central  wavenumber  n-  and  bandwidth  An,  at  pressure 
level  j  with  mean  temperature  Tj  as  in  Equations  (2.21)  and  (2.22). 


B..  = 
ij 

Av.B 

/  V 

J-' 

II 

'dB^ 

\  dB  1 

=  Av, — ^ 

h, 

‘  dt  1 

(2.21) 


(2.22) 
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Channel  Transmittances.  The  channel  transmittances  are  computed  by  applying  the  channel 
spectral  response  a,(v)  to  high-resolution  transmittances  r^<v). 

+  00 

^2.23) 

The  high-resolution  transmittances  were  computed  for  all  wavenumbers  at  each  pressure  level  j 
in  wavenumber  steps  of  An  =  0.01  cm  '  using  FASCOD3  (8].  A  user-defined  atmosphere  was  created 
for  FASCOD3  using  the  pressure  levels  and  mean  temperature/water  vapor  profiles  specified  in  Section 
3  for  the  selected  season/latitude  combination.  The  summer  midlatitude  case  was  chosen  as  the  default. 
FASCOD3  was  run  with  the  seven  major  atmospheric  constituents  including:  H2O,  COj,  O3,  NjO,  CO, 
CH^,  and  Oj.  The  default  FASCOD3  profiles  of  these  constituents  for  the  season/latitude  combination 
were  chosen,  except  for  O2  which  was  artificially  set  low.  The  reduction  of  Oj  level  was  found  to  result 
in  better  behaved  transmittance  functions.  The  details  of  the  FASCOD3  input  are  shown  in  Appendix  A, 
which  contains  the  program  used  to  write  the  input  TAPES  file. 

The  channel  spectral  responses  were  discretized  based  on  their  desired  shape  and  bandwidth  to 
match  the  0.01 -cm  *  step  transmittances.  Figure  1  shows  an  example  for  the  case  of  a  triangular-shaped 
response  function  with  a  full  width  half  maximum  (FWHM)  bandwidth  of  0.05  cm"'. 


Figure  I.  Triangular  spectral  response  array  for  case  with  0.05-cm'^  bandwidth. 
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With  the  spectral  response  of  channel  i  defined  by  normalized  to  a  peak  value  of  one  with  ^^+1 
indices,  the  channel  i  transmittance  for  pressure  level  j  is  found  by  Equation  (2.24).  The  transmittances 
are  centered  at  the  central  wavenumber  v.  of  the  channel  and  summed  with  the  weighting  of  the  spectral 
response  coefficients.  They  are  then  normalized  by  the  sum  of  the  coefficients  to  maintain  the  proper 
amplitude  between  0  and  1. 


I ^j\  ^  +  A v(k  -  JC  /  2)] 

^  a., 

^  ik 


(2.24) 


Computation  of  dt^du.  The  derivative  of  transmittance  of  a  wet  atmosphere  with  respect  to 
precipitable  water  was  computed  by  running  FASCOD3  twice  —  once  for  an  atmosphere  with  water 
vapor  only,  and  a  second  time  with  10%  more  water  vapor.  (In  practice,  FASCOD3  will  not  run  with 
only  one  molecule,  so  a  small  amount  of  Oj  was  added  to  the  atmosphere.)  The  derivative  was  then 
approximated  as  in  Equation  (2.2S)  for  all  of  the  pressure  levels  and  wavenumbers. 


dt  T  —  T 

tv  _  i.iw  V 
du  \AU-u 


(2.25) 


2J  ESTIMATION  OF  RETRIEVAL  OPERATOR  AND  ERROR  COVARIANCE 
2.3.1  Retrieval  Operator 

While  Equation  (2.181  defines  an  ideal  linetu-  forward  radiative  transfer  problem,  a  real  measure¬ 
ment  will  have  a  noise  term  corrupting  the  radiance.  Thus,  we  will  define  the  measured  radiance  vector 
of  the  N  channels  to  be  b^, 

b^  =  Av  +  e  ,  (2.26) 

where  e  is  the  noise  term  assumed  to  have  zero  mean  and  covariance  matrix  N.  The  linear  retrieval 
operator  then  is  defined  to  be  C  which  is  used  to  compute  an  estimate  v  of  the  temperature/water  vapor 
difference  profile  based  on  the  measurement. 

v  =  Cb  (2.27) 

Several  techniques  can  be  used  to  find  a  suitable  retrieval  operator.  For  this  report,  a  minimum 
variance  (MV)  approach  is  used  which  solves  for  C  such  that  the  variance  of  the  estimate  is  minimized. 
Since  the  mean  value  is  assumed  to  be  zero,  this  is  also  a  linear  mean  square  error  (LMSE)  estimate. 
Thus,  we  want  to  minimize  Equation  (2.28)  with  respect  to  C. 
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”1 

=  E\vv'^-vTCb  -b^C^v  +  b^C^Cb  ]  (2.28) 

I  o  o  o  oj 

Taking  the  derivative  with  respect  to  C  and  setting  equal  to  0  yields 
0  =-2£{vi;)  +  2C£{i/j)  . 

Solving  for  C, 

=  £|v(i4v  +  £)^|£|(/lv  +  f)(>li'  +  e)^| 

=  E{vv^A^^ve^}E{Avv^A^  +  Ave^-^sv^A^-^ee^}  ' 

Since  the  noise  term  and  the  solution  difference  vector  v  are  uncorrelated,  and  both  have  zero 
mean,  the  retrieval  operator  becomes  as  in  Equation  (2.29). 

C  =S/1^(aSA^  +  )V)"'  (2.29) 

Here, 

S  =  £{vv^} 

To  avoid  taking  the  inverse  of  a  large  matrix  (the  dimension  of  which  would  be  the  number  of 
channels  which  may  be  in  the  thousands),  the  following  matrix  identity  is  applied. 

SA^(aSA^  +  )V)~' =(a^(V-‘a  +  S"')"‘a^A^"*  (2.30) 

Thus,  the  final  form  for  the  retrieval  operator  is  shown  in  Equation  (2.31). 

C  =  (a^)V“' A  +  S“')"‘ (2.31) 

2.3.2  Error  Covariance 

The  error  in  the  estimate  of  the  temperature/water  vapor  difference  vector  v  is  described  by  the 
covariance  matrix  U,  which  is  computed  as  follows. 
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U  =£■{(!'-  v)(v  -  vf  )  =  £  |(*'  -  Cb^  )(v  -  c  )^  } 

=  E\vv^-Cb  -vb^C^  +  Cb 

I  o  o  00} 

=  s-cas-sa'''c''' +c(asa^  - n)c'^ 

=  S-CAS-SA'^(aSA'^  +  Ny'AS'’' 

+  SA'^(aSA'^  +  Ny\ASA'^  +  N)(ASA'^  +  Ny'  AS'^ 

=  S-CAS  (2.32) 


This  can  be  simplified  through  the  use  of  Equation  (2.31)  and  the  following  matrix  manipulations. 


u 

=  S-(a’'(V-' A  +  S"‘) 

A^N'*  AS 

=  |^/-(a^;V-' A  +  5"' 

=  r(A^N-' A  +  S''r‘ 

(A^iV-'A  +  S"'' 

)-(a^/V-' A  +  S-' 

r‘A^Ar-l/ 

s 

V  j 

'  V 

/ 

=  (A^Af-' A  +  5~')"'[. 

A^/V-'A  +  S-'- 

a^a^~'a]s 

(2.33) 

u 

=  (a^/V'‘ A  +  S"')~‘ 

(2.34) 

The  vector  e  of  retrieval  rms  errors  is  then 

e  =  ^jdiag{U)  (2.35) 

Note  that  throughout  the  analysis,  the  mean  bias  of  the  estimate  is  assumed  to  be  zero;  thus,  the 
expected  error  associated  with  an  instrument  and  given  a  priori  statistics  is  completely  described  by  the 
error  covariance. 
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3.  A  PRIORI  TEMPERATURE  AND  WATER  VAPOR  STATISTICS 


The  a  priori  statistics  used  in  the  retrieval  error  analysis  were  obtained  from  NOAA.  They  were 
derived  from  a  database  of  radiosonde  measurements  [9],  Table  1  shows  the  40  pressure  levels  at  which 
the  values  were  measured.  Temperature  measurements  were  obtained  at  all  40  levels  plus  the  surface, 
while  the  water  vapor  measurements  were  obtained  only  up  to  300  mbar.  Above  that,  the  water  vapor 
concentration  was  assumed  to  be  zero. 


TABLE  1 


Pressure  Levels  Used  in  Analysis  (mbar) 


Pressure 

Level 

Pressure 

Level 

Pressure 

Level 

Pressure  | 

1 

1000 

11 

475 

21 

100 

31 

B 

950 

12 

430 

22 

85 

32 

H 

920 

13 

400 

23 

70 

33 

B 

850 

14 

350 

24 

60 

34 

3.0 

H 

780 

15 

300 

25 

50 

35 

2.0 

B 

700 

16 

250 

26 

30 

36 

1.5 

B 

670 

17 

200 

27 

25 

37 

1.0 

H 

620 

18 

150 

28 

20 

38 

570 

19 

135 

29 

15 

39 

mm 

500 

20 

115 

30 

10 

40 

0.1 

The  statistics  were  obtained  for  six  different  latitude  and  month  combinations,  as  listed  in  Table  2. 
Figures  2  and  3  show  the  mean  temperature  and  water  vapor  profiles,  while  Figures  4  and  5  present  their 
standard  deviations. 
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TABLE  2 


Month  end  Latitude  Statistics 


Month 

Latitude  Region 

January 

0  -  25 

January 

25  -  55  °N 

January 

55  -  90  “N 

August 

0  -  25  “N 

August 

25  -  55  °N 

August 

55  -  90  “N 

16 


21778ft-2 


c 


TEMPERATURE (K) 


Figure  2.  Mean  temperature  profiles  for  January  (top)  and  August  (bottom). 


PRESSURE  (mbar) 


Figure  3.  Mean  water  vapor  profiles  for  January  (top)  and  August  (bottom). 
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PRESSURE  (mbar) 
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4.  EXAMPLE  RETRIEVAL  ERROR  ANALYSES 


4.1  INTRODUCTION 

This  section  presents  two  examples  of  how  the  retrieval  error  analysis  technique  can  be  applied 
to  study  the  expected  retrieval  error  from  infrared  sounding  instruments.  The  first  example  examines  the 
performance  of  the  proposed  Atmospheric  Infrared  Sounder  (AIRS)  instrument  as  an  example  of  the 
technique  applied  to  a  high  spectral  resolution  instrument  with  a  large  number  of  channels.  The  second 
example  shows  the  application  of  the  technique  to  the  GOES-I  sounder,  a  lower  resolution  filter-wheel 
type  of  instrument.  Plots  showing  the  expected  error  are  presented  as  example  results  of  the  analyses. 

A  few  comments  on  the  computational  implementation  are  in  order.  The  algorithm  was  imple¬ 
mented  in  the  IDL  programming  environment  to  take  advantage  of  the  built-in  mathematical  and  matrix 
operations  as  well  as  to  allow  easy  plotting  and  data  visualization.  A  single  parameter  input  file  was  used 
to  define  conditions  for  each  analysis  run.  Also,  an  intermediate  file  was  written  to  disk  containing  the 
A  matrix  for  each  instrument.  This  file  allowed  subsequent  runs  for  the  same  instrument  but  with  different 
noise  levels  to  be  completed  quickly.  Computations  to  create  the  A  mamix  generally  took  1/2  to  1  h  on 
the  Silicon  Graphics  Indigo  R4000,  while  the  N  matrix  formation  and  error  computation  required  only 
15  s. 


4,2  AIRS 

The  proposed  AIRS  is  planned  for  deployment  aboard  NASA’s  Earth  Observing  System  (EOS) 
PM  satellite,  which  is  to  be  launched  in  2(X)0.  The  AIRS  instrument  is  currently  undergoing  design 
changes,  but  for  illustrative  purposes  an  instrument  model  based  on  the  original  specifications  [10]  has 
been  implemented  using  the  error  analysis  methodology. 

The  analysis  was  performed  using  the  August  midlatitude  a  priori  statistics  and  an  A  matrix  formed 
for  the  AIRS  instrument  using  the  information  shown  in  Table  3.  The  specific  channel  locations  were 
taken  from  a  memo  by  NOAA  scientist  David  Wark  [11].  All  channels  were  implemented  except  those 
in  the  ozone  band  (1000-1 100  cm  ')  and  those  above  2400  cm"'.  The  channel  bandwidths  for  the  AIRS 
instrument  are  specified  to  be  «/1200,  where  n  is  the  wavenumber.  The  error  analysis  software,  however, 
requires  the  bandwidths  be  constant  for  a  given  spectral  region.  To  avoid  defining  1229  regions,  the  total 
spectral  bandwidth  was  divided  into  ten  regions,  as  shown  in  Table  4,  with  each  channel  within  a  region 
assigned  the  average  bandwidth  for  that  region. 
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TABLE  3 


Model  AIRS  Instrument  Specifications 


Number  of  channels 

1229 

Channel  set 

WARK 

Spectral  response  shape 

Rectangle 

Number  of  spectral  regions 

10 

TABLE  4 


Model  AIRS  Spectral  Regions  (cm*^) 


Spectral  Region 

Lowest 

Highest 

Bandwidth 

1 

650 

750 

0.58 

2 

750 

875 

0.68 

3 

875 

1000 

0.78 

4 

1100 

1300 

1.00 

5 

1300 

1500 

1.00 

6 

1500 

1700 

1.33 

7 

1700 

1900 

1.50 

8 

1900 

2100 

1.67 

9 

2100 

2300 

1.83 

10 

2300 

2400 

1.96 

The  noise  matrix  N  was  formed  as  a  diagonal  matrix  with  each  entry  being  the  variance  of  the 
noise  equivalent  radiance  (NEN)  for  that  channel.  The  standard  deviation  of  the  NEN  was  calculated 
assuming  an  NEAT  =  0.25  K  for  a  unit  emissivity  blackbody  at  250  K  using  Equation  (4. 1 ). 


‘^NEN 


Av(0.25) 


dT 


r=250K 


(4.1) 


where  Av  is  the  appropriate  channel  bandwidth  and  dBydT  is  the  derivative  of  the  Planck  spectral 
radiance  function,  as  shown  in  Equation  (2.20). 

After  selecting  the  August  midlatitude  case  for  the  a  priori  covariance  matrix  S  and  forming  the 
required  instrument-based  matrices  A  and  N,  the  computation  of  the  error  covariance  matrix  was  per- 
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WATER  VAPOR  MIXING  RATIO  ERROR  (g/kg) 


Figure  6.  AIRS  retrieval  error  for  temperature  (top)  and  water  vapor  (bottom). 
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43  GOES'!  SOUNDER 


The  GOES-I  satellite  [12]  to  be  launched  in  1994  will  carry  an  infrared  sounder  derived  from  the 
HIRS  instrument  on  the  NOAA  polar  satellite.  This  18-channel  filter-wheel  instrument  is  expected  to 
provide  data  that  will  result  in  soundings  comparable  to  those  from  the  polar  satellite,  but  which  will  be 
available  much  more  often  because  of  the  geostationary  orbit. 

Table  5  shows  the  general  instrument  model  as  implemented  for  this  example,  while  Table  6  presents 
the  specifications  of  the  individual  channels.  The  central  wavenumbers  and  bandwidths  have  units  of  cm  ', 
while  the  noise  equivalent  spectral  radiance  (NEAN)  is  given  in  units  of  mW/(m^-Sr-cm  '). 

Figure  7  shows  the  resulting  errors  for  the  temperature  and  for  the  water  vapor  mixing  ratio  when 
using  the  GOES-I  sounder  model.  These  errors  are  significantly  higher  than  those  estimated  for  the  high 
spectral  resolution  sounder  AIRS.  This  type  of  relative  comparison  of  predicted  errors  is  a  useful  appli¬ 
cation  of  this  error  analysis  technique. 


TABLE  5 

Model  GOES'I  Sounder  Specifications 


Number  of  channels 

18 

Spectral  response  shape 

Rectangle 
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TABLE  6 

Model  GOES-I  Sounder  Channels 


Channel  Number 

Central  Wavenumber 

Bandwidth 

NEAN 

1 

680 

13 

0.66 

2 

696 

13 

0.58 

3 

711 

13 

O.wl** 

4 

733 

16 

0.45 

5 

748 

16 

0.44 

6 

790 

30 

0.25 

7 

832 

50 

0.16 

8 

907 

50 

0.16 

9 

1030 

25 

0.33 

10 

1345 

55 

0.16 

11 

1425 

80 

0.12 

12 

1535 

60 

0.15 

13 

2188 

23 

0.013 

14 

2210 

23 

0.013 

15 

2245 

23 

0.013 

16 

2420 

40 

0.008 

17 

2513 

40 

0.0082 

18 

2671 

100 

0.0036 

25 


ESSUR 


TEMPERATURE  ERROR  (‘C) 


Figure  7.  GOES-I  sounder  retrieval  error  for  temperature  (top)  and  water  vapor  (bottom). 


5.  SUMMARY 


The  derivation  and  implementation  of  a  retrieval  error  analysis  scheme  for  passive  infrared 
atmospheric  sounding  instruments  has  been  described.  This  technique  is  based  on  a  linearization  of  the 
radiative  transfer  equation  with  terms  for  temperature  and  water  vapor  explicitly  written.  After  construc¬ 
tion  of  an  instrument  matrix  containing  factors  for  Planck  radiance,  numerical  quadrature  weights,  and 
atmospheric  layer  transmittances,  the  estimated  retrieval  error  is  computed  in  one  matrix  operation. 
Subsequent  analyses  for  various  instrument  noise  levels  can  be  computed  quickly  with  simple  replace¬ 
ment  of  the  entries  in  the  noise  covariance  matrix. 

Examples  were  presented  of  the  technique  applied  to  a  high  spectral  resolution  sounder  with  a 
large  number  of  channels  and  to  a  low  resolution  filter-wheel  sounder.  Predicted  retrieval  rms  errors  of 
temperature  and  water  vapor  for  each  of  these  instruments  were  presented.  The  estimated  error  for  the 
filter-wheel  sounder  was  significantly  higher  than  for  the  high  resolution  sounder  showing  a  relative 
comparison  application  of  the  error  analysis  technique. 

The  retrieval  error  analysis  technique  presented  in  this  report  has  formed  the  basis  for  several 
system  studies  conducted  at  Lincoln  Laboratory,  and  it  has  been  found  to  be  a  useful  tool  in  the  relative 
comparison  of  various  proposed  sounding  instruments  and  design  options. 


APPENDIX  A:  FASCOD3  INPUT 


The  following  FORTRAN  program  was  used  to  generate  the  FASCOD3  input  TAPES  file  which 
was  then  used  to  generate  the  high  spectral  resolution  layer  transmittances  used  in  the  retrieval  error 
analysis.  The  default  conditions  were  a  midlatitude  summer  atmosphere  with  the  seven  major  atmospheric 
constituents  present.  The  input  TAPES  file  was  modified  slightly  to  produce  the  transmittances  for  the 
water-vapor-only  situation  by  zeroing  the  levels  of  all  constituents  except  the  water  vapor.  In  this  case, 
the  level  of  oxygen  was  set  artificially  small  as  it  was  found  that  FASC003  would  not  execute  with  only 
one  molecule  present. 

The  program  below  uses  two  input  files:  TWCM  which  contains  the  meai  .md  covariance  sta¬ 
tistics  of  the  temperature  and  water  vapor  profiles  and  SOUNDPl  which  listed  me  40  pressure  levels 
corresponding  to  the  profile  entries  in  TWCM. 

The  resulting  output  of  FASCOD3  consists  of  the  transmittance  for  all  wavenumbers  for  each  of 
the  layers  in  separate  files.  These  layer  transmittances  then  were  combined  to  produce  the  transmittance 
from  each  layer  to  the  top  of  the  atmosphere,  as  was  needed  for  the  error  analysis  software. 


c 

SOUNDFAS 

c 

c 

J.  Kerekes 

July  22,  1992 

c 

c 

APPLICATION! 

c 

c 

Program  to  write  a  TAPES  input  file  for 

FA5CODE 

c 

using  a  defined  atmospheric  profile  with 

FSCATM. 

c 

Set  up  for  generating  0.01  cm'-l  steps  from 

c 

user  prompted 

range. 

c 

MODEL  2  Midlatitude  summer  defaults 

c 

c 

August  25-55 

N  data  set . 

c 

c 

INPUT  FILES: 

c 

c 

SOUNDPl 

vector  of  pressure  levels 

c 

TWCM 

ten^rature  and  water  vapor 

means  and 

covariance 

c 

c 

matrices 

c 

c 

OUTPUT  FILES: 

c 

c 

TAPES 

input  file  for  FASCOD3  to 

generate  layer  by 

c 

layer  atmospheric  transmittances  over 

desired 

c 

c 

spectral  range 

29 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


+ 

+ 

+ 


LOCAL  VARIABLES: 


aiim 

dvpar 

gconst 

iseaslat 


ibuf 

itempsurf 

maxlevel 

maxmol 

midstart 

numlevel 

nunniol 

niimtenip 

nutnv 

numwv 

outstart 

pm 

q 

rconst 

scov 

startnu 

stopnu 

tbar 

tstar 

vbeu: 

zdiff 


mass  of  air 

wavenumber  step  size  for  interpolated  outpHit 

gravity  acceleration  constant 

index  selecting  season  latitude  con^ination 

of  TWCM  file 

=  1  January  25  -  55  N 

=  2  Jaunuary  55  -  90  N 

=  3  August  55  -  90  N 

=  4  August  25  -  55  N 

=  5  Januau:y  0  -  25  N 

=  6  August  0  -  25  N 

integer  buffer  used  in  reading  TVCM 

index  in  vbar  of  surface  temperature 

maucimum  number  of  atmospheric  levels 

mauciroum  number  of  molecules  in  atmosphere 

initial  unit  number  of  interim  file  sequence 

number  of  atmospheric  levels  used 

number  of  molecules  used  in  atmosphere 

number  of  atmospheric  temperature  layers 

total  number  of  indices  in  vbar 

number  of  atmospheric  water  vapor  layers 

initial  unit  number  for  output  tau's 

aurray  of  pressure  levels 

water  vapor  mixing  ratio 

constamt  used  in  hypsometric  equation 

tenp/wv  covauriance  matrices  from  TWCM 

starting  (lowest)  wavenumber  for  spectral  range 

stopping  (highest)  wavenumber  for  spectral  range 

mean  tenperature  profile  for  selected  iseaslat 

tbau:  adjusted  by  mixing  ratio 

temp/wv  meam  vectors  from  TWCM 

altitude  difference  between  atmospheric  layers 


Other  variatbles  defined  in  FASC003P  User's  Guide 


program  soundfas 

paurameter  ( maxinol=32 ,  maxlevel=67  ,midstart=20 ,  outstart=60 , 
numlevel-4 1 , iseaslat=4 , 

numtenp=40 ,  nuniwv=15 ,  numv=56 ,  itenpsurf =4 1 ,  nuninol=7 , 
dvpar=0.01 ) 
character* 80  cxid 
chau:acter*24  hmod 
chau:acter*2  cmrg 

chairacter*!  jcharp,  jchart,  jchar  (maxlevel, maxmol ) 

real  zbnd(maixlevel)  ,zm(mauclevel)  ,pm(maxlevel)  ,tm(maxlevel) 

real  vmol( maxlevel, mauanol) 

lnteger*4  ibuf (3200) 

real*4  vbau:(numv,6) ,scov(nuinv,numv,6) 
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c 

c 

c 


Prompt  for  input  frequency  range 


write ( 6 ,*)' Enter  starting  wavenumber  in  cm*-l‘ 
read( 5 , * ) startnu 

write(6,* ) 'Enter  ending  wavenumber  in  cm*-l' 
read ( 5 , * ) stopnu 

write(6, ' (a,2f8.2) ' ) '  Range  is  ' ,st£utnu, stopnu 
c 

c  Open  temperatui e/water  vapor  statistics  file 

c  and  read  in  arrays 

c 

open(unit=ll,file='TWCM’ ,status='unknown' ) 
rewind( 11 ) 
do  30  i=l,6 

read(ll,*)  ( ibuf ( j ) , j=l ,3200) 
do  10  j=l,56 

vbao'  ( j ,  i ) =ibuf  ( j )  *  le-6 
10  continue 

incr=56 
do  20  j=l,56 
do  20  )c=l,56 
incr=incr+l 

scov( j ,k, i )=ibuf ( incr) *le-6 
20  continue 

30  continue 

close (11) 
c 

c  Open  and  read  in  pressure  profile 

c 

open  ( unit=12 ,  f  ile= '  SOUNDPl ' ,  status* '  uiiknown ' ) 
rewind ( 12) 
do  15  l=l,numlevel 
read(12,*)pm(l) 

15  continue 

close ( 12) 
c 

c  Open  TAPE5  file 

c 

open ( unit* 1 0 , f ile* ' TAPE5 ' , status* ' unknown ' ) 
rewind ( 1?' 
c 

c  Prompt  for  user  identification  label 

c  Record  1 . 1 

c 

write ( 6 , * ) ' Please  enter  user  label  ( starting  with  $ ) ' 
read(5, ' (a80) ' )cxid 
write(10, ' (a80) ' )cxid 
c 

c  Record  1 . 2 
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ihirac'l 

ilblf4*2 

icntnin»l 

iaersl=0 

i«iiit=l 

iscan=2 

ifiltr“0 

iplot=l 

itest«0 

iatin=l 

cmrg='01  ‘ 

ilas=0 

iins=0 

ixsect=0 

irad=0 

inpts=-l 

npts=-l 

write  ( 10, 102  )ihirac,  ilblf4,  icntiun.iaersl,  iemit,  iscan,  ifiltr 
+  , iplot , itest , latm,  anrg , lias , ins , ixsect , irad , 

+  npts , npts 

format ( 10 ( 4x, il ) , 3x, a2 , 3 ( 4x, il ) , il , i4 , lx, i4 ) 

Record  1.2,1  is  omitted  since  iins=0 


Record  1 . 3 

vl=8tartnu 
v2=stopnu 
sanple=4 
dv8et=0 . 0 
al£ao=0 . 08 
avinass=36 . 0 
dptiiiin=0 . 0 
dptfac-0.001 

write (10,103)vl,v2, sample , dvset , al fao , avmass , dptmin , dptfac 

103  £ormat(8£10.3) 
c 

c  Record  1.4 

c 

tbound^vbar ( itempsurf , iseaslat ) 

sreinisl=0.95 

sremis2-0 . 0 

srenu.s3=0.0 

srre£ll=l .O-sremisl 

srre£12=0.0 

8rre£13=0.0 

write(  10,104  )tbound,srenu.sl,sreinis2,sremis3,srre£ll,srre£12 
+  ,srre£13  • 

104  £onnat(7£10.3) 
c 

c  Skip  records  2  (they  are  not  used  here) 


102 

c 

c 

c 

c 

c 

c 
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Record  3.1  (modified  for  user  defined  profile) 


niodel=0 

itype*2 

iboiax^numleve  1 

nozero=0 

nopmt=l 

ninol=nuniiiol 

ipunch=l 

re=0.0 

hspace=100.0 
f reqbar=0 . 0 
co2mx=0 . 0 

write ( 10, 301 )model,itype, ibmax, nozero, nopmt, nmol, ipunch, 
re , hspace , f  reqbar , co2mx 
format(7i5,5x,4f 10.3) 

Record  3 . 2 

hl=100.0 
h2=0.0 
angle=180.0 
range=0 . 0 
beta=0 . 0 
len=0 

write( 10, 302)hl,h2, angle, range, beta, len 
format ( 5f 10 . 3 , is ) 

Record  3 . 3b 

zbnd(l)  is  the  height  of  level  1 

Level  1  is  at  surface,  Level  41  is  top  of  atnosphere  (0.01  .nbar) 
Use  Hypsometric  equation  from  General  Meteorology  by  Byers 

rconst=8.3143e7 

ainn=28.9 

gconst=980.665 

For  levels  1  through  15  use  water  vapor  mixing  ratio 

zbnd ( 1 ) =0 . 0 
do  23  1=2,15 

tbar=vbar ( itenpsurf-1 , iseaslat ) 
q=0 . 00 1 * vbar ( numv-l+l , iseaslat ) 
tstar=tbar / ( 1 . 0-0 . 6*q) 

zdiff=(le-5)*( (rconst*tstar)/(airm*gconst) )* 

( alog  ( pni(  1-1 ) ) -alog  ( pin(  1 ) ) ) 
zbnd( l)=zbnd(l-l )+zdiff 
continue 


For  levels  16  through  41  eissume  water  vapor  is  zero 


c 
c 

do  25  l=16,numlevel-l 

tstar=vbar ( itempsurf-l , iseaslat ) 

zdif f = ( le-5 ) • ( ( rconst*tstar ) / ( airm*gconst ) ) * 

+  (alog(pra(l-l))-alog(pm(l))) 

zbnd(  1 ) =zbnd(  1- 1 ) -t-zdiff 
25  continue 

c 

c  Do  level  41  using  teinperature  of  level  40  plus  t39-t40 

c 

tst£Lr=vbar  ( 1 ,  iseaslat )  -  ( vbar  ( 2 ,  iseaslat )  -vbar  ( 1 ,  iseaslat ) ) 
zdiff*( le-5)*( (rconst*tstar)/(ainn»gconst) )• 

+  (alog(pm(40)  )-alog(pnv{41) )) 

zbnd(41  )=zbnd(40]-fzdiff 
write(10,303)(zbnd(l),l=l,6) 
write( 10 , 303 ) ( zbnd( 1 ) , 1=9 , 16 ) 
write ( 10,303) (zbnd(l) , 1=17,24) 
write ( 10 , 303 ) ( zbnd( 1 ) , 1=25 , 32 ) 
write(10,303)(zbnd(l),l=33,40) 
write ( 10 , 303 ) ( zbnd ( 1 ) , 1=4 1 , numlevel ) 

303  £ormat(8fl0.3) 
c 

c  Record  3.4,  User  comnents 

c 

inniax=ib(nax 

hmod='  Sounder  Study  4  Sum  Mid' 
write( 10, 304 ) inmax, hmod 

304  format(i5,a24) 
c 

c  Record  3.5,  3.6  repeat  for  each  layer 

c 

c  Get  water  vapor  in  mixing  ratio  ( jchar(*,l)=C) 

c  H20  is  molecule  number  1 

c 

do  33  l=l,numwv 

vmol(l, l)=vbar(numv-l+l, iseaslat) 

33  continue 

do  32  l=nurawv+l, numlevel 
vinol  ( 1 , 1 )  =0 . 0 
32  continue 

c 

c  Get  altitudes  wd  ten^ratures 

c 

do  34  1=1, inmax-1 

zm(l)=zbnd(l) 

tm  ( 1 ) =vbar ( itempsurf-l , iseaslat ) 

34  continue 
zm(41)=zbnd(41) 

tm( 4 1 ) =vbar ( 1 , iseaslat ) - ( vbar ( 2 , iseaslat ) -vbar ( 1 , iseaslat ) ) 

jcharp='A' 

jchart='A' 
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c 

c 

Set  up  models  to  use 

for 

other 

molecules 

c 

»•****•*  »******niodel=6 

is 

US  Standeurd 

c 

****»*»******»*jn£xiel=2 

is 

summer 

midlatitude 

c 

do  37  l=l,numlevel 

jchar(l,l)='C' 
do  35  k=2 ,  nmol 
jchar ( 1 , k ) = ' 2 ■ 
vmol ( 1 , k ) =0 . 0 

35  continue 
c 

C  MODIFICATION  FOR  02  HERE 

c  (reduced  02  level  led  to  smoother  tau's) 

c 

jchcur(l,7)  =  'C' 
vmol(l,7)=0.001 
37  continue 

do  36  1=1, inmax 

write(  10, 305)zm(l)  ,pni(l)  ,tm(l) , 

+  jcharp,  jchart,  (  jchar(l,k)  ,k=l,nniol) 

write ( 10, 306 ) (vmol (l,k) ,k=l,  nmol) 

36  continue 

305  £onnat(3£10.3,5x,2al,3x,28al) 

306  £onnat(8£10.3) 
c 

c  Record  9.1  (repeat  for  all  layers,  end  with  negative  dv) 

c 

dv=dvpea' 

vl=vl 

v2=v2 

jatiit=0 

i4pt=0 

iunit=10 

n£ils=l 

npts=-l 

do  91  i£ilst=l,ininax-l 
junit=midstart+i£ilst-l 

write ( 10,901 )dv,vl,v2, jemit,i4pt,iunit,i£ilst,n£ils, 

+  junit,npts 

91  continue 

c 

c  end  with  negative  dv 

c 

dv=(-1.0)*dv 

write ( 10 , 901 ) dv , vl , v2 , jemit , i4pt , iunit , i£ilst , nif ils , 

+  junit,npts 

901  £ormat(3£10.3,2i5, 15x,5i5) 

c 

c  Record  11.1 

c 
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write( 10, ' (a30) ' )cxid 


c  Record  11.2a  (repeat  11.2a  emd  11.3a  for  number  of  layers, 

c  end  with  negative  vl) 

c 

c  Records  11.2a  and  11.3a  are  used  to  produce  ascii  files  o 

c  the  trcuismittances .  Other  "plotting"  vari^d5les  aire  arbitrary, 

c 

xsi2e=9 . 0 
delv=5.0 
nuinsbx=5 
noendx=l 
lskipf=0 
scale=l . 0 
iopt=0 
i4p=0 
ixdec=0 
c 

c  Record  11.3a 

c 

ymin=0 . 0 

yinax=0.9 

ysi2e=6.5 

dely=0 . 1 

numsby=5 

noendy=0 

idec=2 

jemit=0 

jplot®0 

logplt=0 

jhdr=l 

jduniny=0 

jout=3 

do  12  l=l,ininax-l 
lfile=midstart+l-l 

write(10, 112)vl,v2,xsi2e,delv,nunisbx,noendx,lfile, 

+  Iskipf, scale, iopt,i4p,ixdec 

jpltf l=outstart+l- 1 

write(  10, 113)ymin,ynuuc,ysi2e,dely,nuinsby,noendy,idec,  jemit, 

+  jplot,logplt, jhdr, jdumny, jout, jpltfl 

12  continue 

112  fonnat(4f 10.4,4i5, f 10.3,i2,i3,i5} 

113  forTnat{2fl0.4,2fl0.3,6i5,i2,i3,i2,i3) 
c 

c  End  with  negative  vl 

c 

vl=(-1.0)*vl 

write ( 1 0 , 1 12 ) vl , v2 , xsize , delv, numsbx, noendx, If ile , 
t  Iskipf , scale, iopt,i4p,ixdec 
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End  record 


c 
c 

write( 10, ' (a5) ' ) '%  End' 

close(lO) 

stop 

end 
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LIST  OF  SYMBOLS 


A  Radiative  transfer  matrix 

C  Retrieval  operator  matrix 

N  Sensor  noise  covariance  matrix 

S  A  priori  covariance  matrix  of  temperatures  and  water  vapor  mixing  ratios 
U  Retrieval  error  covariance  matrix 

£  {•)  Expectation  operator 

B  Radiance  emitted  by  layer  of  atmosphere  in  a  particular  spectral  channel 
Spectral  radiance  emitted  by  layer  of  atmosphere 
Cj  Planck  function  constant  =  1.191062  x  10'^  mW/m^-Sr-fcm  ’)"* 

Cj  Planck  function  constant  =  1.438786  K/cm"' 

R  Upwelling  radiance  received  by  satellite  in  a  particular  spectral  channel 
T,t  Temperature  of  layer  of  atmosphere  (Kelvin) 

X  Layer  in  model  atmosphere 

b  Vector  of  radiance  differences  between  mean  and  true  value 

bg  Vector  of  radiance  differences  between  mean  and  measured  value 

e  Vector  of  retrieval  error  standard  deviations 

V  Solution  vector  of  temperature  and  water  vapor  differences  between  mean  and  true  value 

V  Estimated  vector  of  temperature  and  water  vapor  differences  between  mean  and  measured  value 

£  Vector  of  sensor  noise  radiance 

a.(v)  Channel  i  spectral  response  function 

Channel  i  spectral  response  at  index  k 

p  Atmospheric  pressure  (mbar) 

q  Water  vapor  mixing  ratio  (g/kg) 

u  Precipitable  water  (cm) 

V  Wavenumber  (cm'*) 

T  Transmittance  from  a  particular  layer  to  the  top  of  the  atmosphere 

Tp  Transmittance  to  top  of  atmosphere  with  no  water  vapor  present 
Transmittance  to  top  of  atmosphere  with  only  water  vapor  present 
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